Mesoscopic theory for size- and charge- asymmetric ionic 
systems. I. Case of extreme asymmetry 

A.Ciach 1 , W.T. Gozdz 1 and G.Stell 2 

1 Institute of Physical Chemistry, Polish Academy of Sciences 01-224 Warszawa, Poland 
2 Department of Chemistry, State University of 
New York Stony Brook, NY 11 794-3400, USA 
(Dated: February 2, 2008) 

Abstract 

A mesoscopic theory for the primitive model of ionic systems is developed for arbitrary size, 
A = cr+/cr_, and charge, Z = e+/|e_|, asymmetry. Our theory is an extension of the theory we 
developed earlier for the restricted primitive model. The case of extreme asymmetries A — > oo 
and Z — > oo is studied in some detail in a mean-field approximation. The phase diagram and 
correlation functions are obtained in the asymptotic regime A —* oo and Z — » oo, and for infinite 
dilution of the larger ions (volume fraction n p ~ \jZ or less). We find a coexistence between a 
very dilute 'gas' phase and a crystalline phase in which the macroions form a bcc structure with 
the lattice constant 3.6cr + . Such coexistence was observed experimentally in deionized aqueous 
solutions of highly charged colloidal particles. 
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I. INTRODUCTION 



For many years, theoretical studies of phase behavior in ionic solutions have been focused 
mainly on the special case of the restricted primitive model (RPM), in which half of equal- 
sized charged hard spheres carry positive charge and half carry negative charge of equal 
magnitude, with the ions assumed to be dissolved in a structureless solvent jll Q, S, Q|. 
Even the simplest real ionic solutions have some degree of size asymmetry, but with some 
notable exceptions 0,0,121 a common tacit assumption has been that the effects of weak and 



moderate asymmetry in both size and charge is not important to phase behavior. Recently, 



the size and c harg e asymmetric case has drawn increasing attention 
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Most extensions beyond the RPM are based either on the Debye-Hiickel 



theory and Poisson-Boltzmann equation, or on the mean spherical approximation. These 



theories 
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as well as simulations 



0, S a, 12, 0, 



are typically 



limited to the case of small differences in sizes and charges. Only in very recent simulations 



have moderate 



121 ] and large 



00 



asymmetries been studied. Moreover, these theories 



are all "classical" [i.e., mean- field-like] , and none of them are designed to describe the special 
Ising-like behavior that is known to characterize the primitive model in its critical region. 
The development of a theory that does describe that behavior was sketched by one of us 



in Ref 



20 and further developed in Refill] and Refl^, A field-theoretic method that also 



yields the correct Ising-like behavior was given by Ciach and Stell in Ref. |23j and further 



developed in our subsequent papers 



000 



l27| . The mesoscopic theory described in 



this paper is an extension of the theory given there, and reduces to it when applied to the 
restricted primitive model. 

In the case of extreme asymmetry (charge- and diameter ratios between the two kinds of 
ions tend to infinity) the PM potentials describe highly charged colloidal particles suspended 
in a structureless solvent containing one kind of counterions and no coions. The physical 
properties of such a system are significantly different than those of the usual electrolytes. 
Highly asymmetric systems exhibit an interesting phase behavior which is neither fully 
described nor understood, but it is clearly quite different than that of the RPM. In particular, 
ormation of a colloidal bec crystal with large inter-particle separation coexis ting with voids 
28 . 29I I30L v arious crystals formed by oppositely charged colloidal particles |3l| . and other 



anomalies [sol . 



HQ 



have been observed. The experimental findings suggest the existence of 
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effective attractions between like-charged macroions as a possible explanation of the observed 



phase behavior. The classic Derjaguin-Landau-Verwey-Overbeek (DLVO) theory 
however, predicts purely repulsive interactions between the like-charged colloids. In recen ; 
approaches geometrical effects, as well as fluctuations and correlations are included H36l . 137 . 

4311 . The approaches are based e.g. on integral equations [38|, 39], 

H 0, Q, and variational methods y. The effective attractions 
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38, 



39, 



density functionals 
appear in a modified DLVO theory 



44( | . and will also result from the explicit inclusion of 



various effects, such as 'charge regulation' |38J], excluded volume ('Coulomb depletion') |4l| . 
metastable states 3] etc. On the other hand, 'volume terms' considered in Ref. Q lead 
to a phase separation for purely repulsive interactions. Some of the mentioned works are 
questioned by authors of the other papers. It goes beyond the scope of our paper to discuss 
the above approaches in more detail; extensive lists of recent works and discussions can 



be found in Ref 



37 
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46l Despite impressive progress, the experimentally observed 



void - bcc crystal coexistence has not been predicted, and the issue is still controversial. In 
theoretical approaches to colloidal systems one typically assumes extreme size asymmetry 
between the macro- and microions, and the methods differ from those developed for the 
RPM or for the PM with a small asymmetry. 

In principle it should be possible to analyze the evolution of phase diagrams when the 
size- and charge ratios increase from unity to infinity. To achieve this goal one needs a 
theory applicable to arbitrary size- and charge asymmetries for the PM potentials. Within 
the context of Ornstein-Zernike formalism, one can go quite far in obtaining the general 
structure of such a theory, from which a number of important general results follow, such 
as the relation between the charge-charge and density-density correlation lengths, which 



shows that they must diverge together in the asymmetric case [7J]. However, quantitative 
results for the thermodynamics and structure of systems of asymmetric ions are very sensi- 



tive to approximations and assumptions 



In fact one often needs to know 



the results to make proper assumptions, i.e. to identify the physical effects that have to 



be espncitly taken into acconnt (association Q, ^order Z one- M andclnster U forma 



n 



tion, 'charge regulation' |38j, 'Coulomb depletion' |4l|, 'volume terms' |45( etc.). Recently 
dev elop ed field theory for asymmetric ions, based on the Hubbard- Stratonovich transform 



47 



42 1, is elegant and in principle exact. In practice, however, the phase eqilibria and 



correlation functions can be obtained by using different approximate methods. In Ref 



47 
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the size of ions is taken into account through the single cutoff in the Fourier space. Since 
in ionic systems the dominant fluctuations are short-range charge-density waves [5^], the 
short-distance properties of the system are important, and this approximation may lead to 
inaccurate results. Moreover, the effect of size asymmetry cannot be studied in the theory 
with a single microscopic length. In Ref. |49| the hard spheres are taken into account more 
directly, and the charges are smeared inside the spheres to regulate the Coulomb potential. 
Formal expressions and relations are derived for arbitrary asymmetry, but in the general case 
they are very complex (and depend on the smearing function) and the author restricts the 
analysis of the thermodynamics and structure to the case of equally sized ions. Reliable the- 
ory allowing for a determination, with a reasonable effort, of phase equilibria and structure 
in the case of arbitrary asymmetry between the ions has not been developed yet. Therefore 
the crossover between the case of full symmetry and the case of extreme asymmetry is an 
essentially unexplored problem. The PM in the crossover region might be an appropriate 
model for ionic liquids, and it certainly deserves attention. 

In a tractable theory simplifying assumptions and approximations are unavoidable. The 
key issue is to identify the degrees of freedom relevant for phase transitions and critical phe- 
nomena (i.e. along the spinodal lines), and to develop a theory which takes them into account 
correctly, with the irrelevant degrees of freedom treated in an approximate way. In order 
to describe phase transitions where ordering occurs at the length scales large compared to 
molecular sizes, a coarse-graining procedure, leading to the Landau-Ginzgurg- Wilson (LGW) 
approach, has been introduced. The basic assumption of the LGW theory is that for macro- 
scopic phase separation the short-wavelength fluctuations, and hence the precise form of 
correlations at distances r pa a are irrelavant. In the case of simple fluids a correlation 
function can be thought of as being the sum of two pieces- the piece that is on the scale of 
the distance between particles plus the piece that is on the scale of the correlation length, 
which is arbitrarily large close to the spinodal, and it is only the latter piece that deter- 
mines universality class and critical exponents. One can neglect the short-range behavior of 
correlation as long as one is in a critical region, [but only then]. 

In the coarse-grained description one considers deviations from random distributions of 
molecules, and it is important to include the dominant, most probable fluctuations. In simple 
fluids these correspond to macroscopic separation, i.e. to fluctuations with the wavenumber 
k — > in Fourier representation. Because like-charge ions repell, and oppositely-charge ions 
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attract each other, in ionic systems charge-ordered clusters, where positive- and negative 



charge ions are the nearest-neighbors, are observed in real space 



23 



12 



5lj . In Fourier represen- 
. Thus, the fluctuations 



tation the dominant fluctuations are charge- density waves 
associated with charge ordering in periodic structures should be included in the coarse- 
grained description. 

The idea of coarse-graining was sucessfuly extended by Brazovskii [5^ and others to soft- 
matter systems (liquid crystals, microemulsions, diblock copolymers), where microphase 
separation occurs, i.e. periodic phases with a mesoscopic period of density oscillations may 
become stable. A mesoscopic period means a period of order of several molecular diameters 
or larger. In this case one expects that the phase equilibria should be qualitatively correctly 
described, provided that the fluctuations on the length scale corresponding to the ordering 
are included. Again, the correlation function consists of a short-distance piece and of the 
piece that oscillates on the mesoscopic scale and decays on the scale of the correlation length, 
which is arbitrarily large close to the spinodal. It is the latter piece that determines the 
phase transitions, as in simple fluids. The separation into the short- and long-distance 
pieces of the correlation function can be conveniently done by a pole analysis in the complex 



Fourier space 



53, 



54. 



55, 



561 ] . The dominant pole (or a pair of complex conjugate poles) 



with the smallest imaginary part determines the asymptotic large-distance behavior. It 
turns out that this dominant pole (or a pair of poles) decribes quite correctly the correlation 
'unction down to the second maximum for short-range and for Coulombic interactions 
[n the coarse-grained mesoscopic theories the remainig poles of the correlation 
functions are neglected. Note, however that down to the second maximum in the correlation 
functions the neglected poles lead to a small correction to the correlation function, and the 
results of the Landau-type, mesoscopic theories work well for such distances, although for 
distances r < a they are meaningless. To conclude, if one is interested in the vicinity of 
the spinodal line and in the large-distance part of the correlation functions [but only then], 
and one wants to take into acount the possiblity of ordering at distances corresponding to 
the second maximum of the correlation function or larger, one can consider the mesoscopic 
Landau-Brazovskii theory. 



Here we propose to extend the mesoscopic field theory introduced for the RPM in Ref 



23 



to the case of arbitrary size- and charge asymmetry. The results of our field theory agree 
with simulations in continuum- space RPM Q,Q> on the sc and the 
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fee lattices 



26 



and on finely discretized versions of the former l26l 127 ] , and also in 



the presence of additional short-range attractive 



63j | and repulsive |26j, |56j interactions. 



Moreover, the electrostatic free energy has a correct behavior for low densities, and the 



exact result in the Debye-Hiickel limit is correctly reproduced [25]. So far no example of 
qualitatively wrong predictions of the mesoscopic theory has been found, although in some 
cases (including the RPM in continuum space) the effect of fluctuations has to be properly 



taken into account 



26 



27 



561 ] : this can be done systematically in the perturbation theory. 



Foundations of the mesoscopic theory do not depend on the symmetry properties between 
the ionic species. On the basis of the results obtained for different extensions of the RPM 
one can hope that an extension to the case of arbitrary asymmetry will also result in a 
predictive theory yielding correct results on a semiquantitative level. 

In this work we introduce the general framework of the mesoscopic theory for the PM 
(sec. 2). Next, in sec. 3, we focus on the case of extreme asymmetry, which turns out to be 
particularly simple. We find the phase behavior and compare it with experimental results 
for highly charged colloidal particles in salt-free water [3, 0, 13 ■ The agreement is very 
good. We also derive the correlation functions for extremely asymmetric case and show 
their forms for various thermodynamic states. We obtain monotonic decay of correlations 
for very dilute system, and results consistent with electric double layer formation for less 
dilute systems. Near the transition to the bee structure the double layer becomes denser and 
thiner. We find a pronounced maximum of the colloid correlation function at distances that 

n 

agree with experimentally observed ordering |3i| • At such distances the clouds of counterions 
around the particles do not overlap. Our results indicate that the theory developed for 



arbitrary asymmetry leads to qua 



- fully symmetric (RPM) 



24 



itatively correct predictions in two opposite limiting cases 



25 



26 



561 ] and extremely asymmetric. Hence, we can 



expect qualitatively correct results in the crossover region as well. The results in the case of 
arbitrary asymmetry will be described elsewhere. 
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II. MESOSCOPIC THEORY FOR THE PM 



A. Coarse-graining procedure 

We consider the PM electrolytes with the diameter and charge ratio between the large 
and small ions 

o+/o-_ = A and e + /|e_| = Z (1) 

respectively (without a loss of generality we assume e_ = — |e_|,e+ = |e+|). For Z — A = 1 
the model reduces to the RPM, and for Z, A — > oo the model describes highly charged colloid 
particles and point-like counterions with a small charge. In the PM the interaction potential 
of a pair a, (5 — ± is infinite for distances smaller than the sum of radii, 

CT a /3 = (<T a + 0/?)/2, (2) 

i.e. we assume hard-core repulsions. The electrostatic potential V a p(r 1 — r 2 ) between the 
pair of ions a, (5 is 

U) = ^-^)> (3) 

where D is the dielectric constant of the solvent (water). The ^-functions above exclude the 
contributions to the electrostatic energy coming from overlapping hard spheres. 

In our field-theoretic, coarse-grained approach, we consider local instantaneous densities 
of the ionic species, p Q (r), i.e. we specify the numbers of ions of both kinds per mesoscopic 
volume dr. For given densities p a (r) precise positions of the ions can be different, and the 
probability density p that the local densities assume a particular form p+(r), p_(r) is given 
by 

P[p a (r)] = S- 1 ! e-W**\ (4) 

J Sp 

where (3 = (ZcT) -1 , and where T and k are temperature and the Boltzmann constant respec- 
tively. By J s we denote an integral over all microscopic states S p compatible with the chosen 
densities p+(r), p_(r), and by E(S P ) we denote the energy of the corresponding microstate. 
The energy of the microstate S p can be written in the form E(S P ) = C/[p + ,p_] + AE P (S P ), 
where ?7[p + ,p_] = J s E(S p )/Af is the mean energy for fixed densities p+(r), p_(r), and 
M — j s is the number of all microscopic states compatible with p+(r), p-(r). We assume 
that for all microscopic states compatible with the given local densities the energy of the 
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whole system is approximately the same, so that j3AE p (S p ) -C 1. Hence, 

e -ms P ) = e -pu[P + , P -] [i _ (3AE P (S P ) + l -{^E p {S p )) 2 + ...] , (5) 
and the probability (J1J) can be written in the form 

p = E- 1 e- pu[p+ ' p - ] [Af + corr] . (6) 

The correction term is proportional to f s (/3AE p ) 2 <C M and will be neglected. Finally, we 
assume that for particular fields p + (r),p_(r) the electrostatic energy is given by 

U[p+,P-} = \ J dr i J dr 2 p a (r 1 )V a / 3 (r 1 - r 2 )p /3 (r 2 ), (7) 

where the summation convention for Greek indeces is used. The fields p a (r) for which 
U — *■ oo occur with the probability p —>■ 0. Hence, in macroscopic regions the charge 
neutrality condition 

J drp + (v)e + = J drp_(r)|e_| (8) 

must be obeyed. One can easily verify that when (jHj) is satisfied for uniform fields p a (r) = 
const, then U[p + , p_] = 0. Due to thermal motion the charge neutrality can be violated in 
mesoscopic regions containing a small number of ions. The energy ((7j) associated with local 
deviations from the charge neutrality remains finite. 

When AE(S P ) can be neglected (i.e. for all microscopic states compatible with p a (r) 
the energy is approximately the same), we can use the Boltzmann formula M = exp(/3TS), 
where by S we denote entropy. In an open system the probability is also proportional to 
the activities exp[/3(p + N + + p_N_)], where N a = J drp a (r) is the number of ions of the 
species a, and the chemical potentials p a are not independent - they have to be consistent 
with the requirement of the charge neutrality (JHJ). The above discussion shows that the local 
instantaneous densities assume the form p + (r),p_(r) with the probability density given by 

p = E- l eM-P^ MF [P + ,P-]), (9) 

where Q MF [p + , p_] is the grand potential in the system where the local concentrations of 
the two ionic species are constrained to be p+(r), p_(r). Next we assume that the entropy 
is determined by the hard-core reference system with the Helmholtz free energy Fh = —TS, 
and fl MF [p + , p_] is assumed to have the form 

fi MF [p+,p_] = F h [p + ,p_] + U\p+,p-] - J drp aPa (v). (10) 



For the reference system we assume the local density approximation Fh[p+,p~-] = 
J drfh(p+(r), p_(r)). The fh consists of the ideal-gas contribution plus the excess free- 
energy density of hard- spheres with different diameters f^ x . For example, the Percus-Yevick 



approximation for hard-sphere mixtures 64| can be adopted. Because of the above assump- 
tion, packing effects of hard spheres cannot be described in our theory, and in the present 
form it is not applicable to very high densities. In principle, extensions beyond the local 
density approximation are also possible. 

In the field theory introduced above the physical quantities are obtained by averaging over 
all fields p+,p_ with the Boltzmann factor (JHJ). The average densities and the correlation 
function are respectively given by 

W- 1 I Tin. I Tin P -W MF [P + ,P-] 



and 



<p Q (r)> = H- 1 J D P+ J Dp-.e-» l " ^ Pa (r) (11) 

G aP (r,r>) = ( Pa {v)p p {v')) - (p a {r))(p p {v>)) (12) 

(p«(r) P/9 (rO> = S- 1 J D P+ J Dp_e-^ MF ^p a (r)p p (r>) } (13) 

Dp + [ Dp_e-^ MF ^~\ (14) 



with 



and 



The grand potential Q is 

-/3fi = logS. (15) 

In practice we are not able to evaluate the functional integrals in Eqs. (fTTj). (JT3~|l and (JHJ), 

and we need to make approximations. In the simplest, mean-field (MF) approximation the 
average values of the local densities are approximated by their most probable values, po a , 
and the grand thermodynamic potential is approximated by the minimum of Q MF [p + , p_] 
at p a = po a . 

As convenient thermodynamic variables we choose dimensionless number density of all 
ionic species, s, and dimensionless temperature T* = where 

* = £(pk+P&-), ^ = 4^' ( 16 ) 

and 

P* = a%_p a . (17) 



Here and below as a length unit we choose . Because of the charge neutrality, 



Pa- = Z P 



(18) 



the volume fraction of ionic species is 



C = 6 {p0+a + + P -"-) = (1 + A)3(1 + Z) S - 



(19) 



Let us study the form of fl MF in more detail. For small deviations Ap* (x) = p„(x) — Po Q of 
the local densities from their most probable values the grand potential (jl(J|) can be expanded 



about its value Qq IF at the minimum, 



aq mf = n MF - n™ F = n* IF + nfj. 



(20) 



Here Q^ F denotes the Gaussian part of the functional. In Fourier representation we have 

1 f dk 



pn 



MF 



(2tt) 



r A&(-k)C^(k)A#(k) 



(21) 



where Ap* (k) is the Fourier transform of Ap* (x), and the wave numbers are in er + i units. 
The second functional derivatives of fl MF , C®Jk), consist of two terms, 



C ap (k) = a a p + pV a p(k). 



(22) 



The first term is given by the corresponding second derivative of (3fh taken at p* = pQ a . 
The second term in C^(k) is the Fourier transform of the potential (p3j). and we find 

e a e/3 in cos(kr a /3) 



PV a p{k) 



e + |e_| k 2 

where r a p = a a p/<j + -. The remaining part of the functional has the expansion 



(23) 



pa 



MF 
int 



dr 



^Ap;(r)Ap;(r)A P ;(r) + ^Ap;(r) Ap^(r) A P ;(r) ApJ(r) 



where 



<9cw 



and 



da, 



d P ; da ^ s dpi ' 

and the derivatives are taken at p* = p^ a . For pair-potentials and for the local-density 
approximation for the reference system, flf^f is strictly local. The explicit forms of the 
coefficients depend on the reference system. 



a/37 



+ ..., (24) 



(25) 
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B. Phase transitions 



Let us focus on the boundary of stability of the uniform phase. The uniform phase is 
unstable with respect to fluctuations p* (k) when the second functional derivative of Q MF 
is not positive definite, i.e. det C^fc) < 0. The temperature at the instability with respect 
to the k-mode is thus given by 

detcyk) = 0. (26) 

Boundary of stability with respect to the deviations Ap* oc cos(r • k) from the densities 
PQ a corresponds to k = k b such that the Eq.(|2*o^l is satisfied first when the temperature is 
decreased. For fixed pQ a the boundary of stability is thus given by the maximum of T*(k) 
obtained from (|26|). therefore k;, can be determined from the equations 

Solutions of the set of equations (J2tj|) and ()27|1 give both, the wave vector of the critical 
fluctuations k&, and the spinodal line in the phase space (s,T*). For temperatures higher 
than at the spinodal line the randomly chosen instantaneous local densities are most probably 
uniform. For lower temperatures, however, the randomly chosen instantaneous densities 
most probably have a form of planar waves with the wave vector k&, or of linear combinations 
of such waves with different orientations of the wave vectors. The amplitudes of the density 
waves of the ionic species and the order of the associated phase transition depend on the 
form of flint- m the next section we shall find the phase transition in the case of extreme 
asymmetry in MF. 



C. Correlation functions 



Let us consider the structure of the disordered phase, i.e. the correlation functions for 
the density deviations from pg a . In our mesoscopic theory, especially in the local density 
approximation, the correlation functions defined in Eg. ([13)1 are meaningful for distances 
larger than a a p, and in principle we can only expect a semiquantitative agreement with 
results of exact theories or simulations for large distances. The pole analysis of the correlation 
functions in Fourier representation shows that in the mesoscopic theory for the RPM only 
the dominant poles, characterizing the long-distance behavior, are present 65j]. On the other 
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hand, the dominant poles y ield quite accurate results down to the second maximum of the 



correlation functions 



54, 



57(, and for such distances we can expect semiquantitatively correct 
results in the colloid limit as well. However, the functional integrals in Eq. (J13|) cannot be 
calculated exactly. In practice we are able to calculate G a p{r) in a perturbation expansion 
in 72n,m- I n the Gaussian approximation, AQ MF = , i.e. with the term fif^f in Eq. (j2(J|) 
neglected, the correlation functions (jl3j) can be easily calculated by inverting the matrix of 
second functional derivatives of Q MF . In Fourier representation we have thus 

G^(k)= [C°(k)]^. (28) 

The functions analogous to pair distribution functions are related to G a p (Eq.(JT3J)) by 

gUr) = Gg>_M& + 1 . (29) 

POaPo/3 Pa 

Beyond the Gaussian approximation we expect corrections to the correlation functions. 
Their relevance in different thermodynamic states will be studied in future works. 

III. CASE OF EXTREME ASYMMETRY 

The above model can be applied to a suspension of highly charged colloid particles in 
a (salt-free) solvent containing one kind of point-like counterions. In the case of extreme 
asymmetry the reference system corresponds to a mixture of hard spheres and point-like 
species, with the densities constrained according to Eq. (jl8|) . For A — ► oo the volume fraction 
C reduces to the volume fraction n p of the large species, and in the asymptotic regime Z — > oo 
(such that Z/X 3 -> 0) we obtain ((see (fT3jl . IflBjl and IflBjl ). 

C = n p = 8s/Z. (30) 

Hence, finite values of the number density, s = 0(Z°), correspond to infinite dilution of 
hard spheres for Z — > oo. At infinite dilution a hard-sphere system can be approximated 
by an ideal gas. The smaller ions are point-like in the limit A — > oo, and also behave as an 
ideal gas. Thus, for A, Z — > oo we can assume that the reference system is just a mixture 
of ideal-gases. For a mixture of ideal-gases a a p = 5^ '/ *Po a , and from Eqs. (fT7)J) and (JTBJ) we 
obtain 

1 71 1 71 

a = = — , a_i__i_ = = — Z, a+_ = 0. (31) 

Pi- 6s' ++ pl + 6s + { } 
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Strictly speaking, in our analysis we assume that A — > oo first, and next we consider the 
asymptotic behavior of Z —>■ oo with s = 0(Z°). The above asymptotic behavior will be 
referred to as the colloid limit. From Eq. (J23|) we easily find that in the colloid limit the 
electrostatic potentials are 



pv ++ (k) = O(Z), (3V + .{k) = 0(Z°), pv_{k) = 0{Z- 1 ). 
Hence, in the colloid limit we obtain 

C ++ (k) = ZC p {k), C+4k) = f3V + 4k) = 0{Z% £__(*) 



f + 0(*- 1 ) 

OS 



and in turn 



where 



detC a p(k) = Z-C p (k) + 0(l), 

OS 



C p (k) 



71 



1 | 24oos(2fc) 
s k 2 



(32) 



(33) 



(34) 



(35) 



A. Stability of the disordered phase 



Instability of the disordered phase, in general given in Eq. 



to 



[C p {k)^ + O(Z- 1 )]=0 1 



for Z ^> 1 is equivalent 



(36) 



and for Z — > oo the latter equation is satisfied when C p = 0. From Eq. (J33|) we obtain that 
if C p = 0, then C ++ (k) = for arbitrarily large Z. This means that in the considered 
asymptotic regime of Z — > oo and n p = 0(Z^ 1 ), the fluctuations Ap + (k) can destabilize the 
uniform phase. The line of instability of the uniform phase with respect to these fluctuations, 
given by (|26|) and (|27|). assumes the form 



C p (k) = = dC p {k)/dk. 



The spinodal line is given by the explicit expression 



t; (s ) = - WM S , tan(2fct) = i 



(37) 



(38) 



and we find kb ~ 1.23 in a^l_ units. 

In order to determine the phase transition associated with the spinodal (|38|). let us 
consider the asymptotic behavior of Q int (Eq. (p0|l) for Z — >• oo with s = 0(Z°) and A — > oo. 
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For the reference system corresponding to a mixture of ideal gases the only nonvanishing 
coefficients in Eq.(J2UJ) are (see Eq.(|3*T|)) 



1 



Pol 



After using Eqs. (fTOj) and (fTHj) we obtain 

Z 



/3AQ 



MF 



Z 2 A 3 
3! 



A P ; 3 (r 



Z 3 A 4 
4! 



Ap;(k)C7 p (A;)Ap;(-k) + 
A P ; 4 (r) + ....l +O(Z ), 



(39) 



(40) 




where 

- 4 ---(s)- ^ 2 (<sJ- < 41 > 

and by 0(Z°) we denote all remaining contributions to AQ MF with the integrands that 
remain finite or tend to zero for Z — > oo. All integrands proportional to Ap*_ turn out to be 
0(Z°). For Z — > oo we neglect such terms compared to those given in Eq. (|40jl . Note that we 
again come to the conclusion that the phase transition in the colloid limit is determined only 
by the macro ion-density fluctuations. In the second step we rescale the field, Ap*^_ = ip/Z, 
and the functional, AQ[Ap + , Ap_] = £l r [if;]/Z, and we obtain 

/».#] = \j m)C,(k)i\-k) (42) 



<->■ Jki </k 2 </k 3 ■ ^ 

4 / / / / *<x»ri^.). 

4! Al Ad A.< A4 



'ki ^k2 «/k3 ^k4 

where we truncated the expansion in the field at the fourth order term. Because the cubic 
term is present, the transition to the phase with periodic ordering of the particles is first 
order in MF. A similar functional was already studied by Leibler in the context of block 
copolymers 0], and we can directly use his results. In order to find the stable structure 
one considers ip of a form of linear superpositions of n planar waves with the wave vectors 
\t 3 b having different orientations (with j = 1, ...,n and \\t J b \ = kb), 

i n 

^ k ) = -fr E [ 5 ( k - k >- + 5 ( k + k M ' ( 43 ) 



3 
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where w n , u>* are complex conjugate and w n u>* = 1. For ^(k) of the form (jjHj) the functional 
Q r (|4*2j) per volume V can be written as 

PQr/V = C p {k b )$ 2 - a 3 $ 3 + a 4 $ 4 , (44) 



66 for several 



where the geometric factors atj depend on n and have been found in Ref. 
structures. For metastable structures fl r assumes local minima, 

d[3n r /d<S> = 0, (45) 

and the stable structure corresponds to the global minimum. At the coexistence of the 
disordered phase with the periodic structure 

(3Q r = 0, (46) 

and the actual phase transition occurs when the above equation is satisfied for this phase 
which corresponds to the global minimum of (3Q r ; another words, for the phase which 
becomes stable, Eq. (|46jl is satisfied at the highest temperature. From (f4*4*j) . (I43J) and (|4*E|) we 
obtain the transition line as C p {k b ) = a\/{Aa^). From the results of Refl66l it follows that 
in MF the disordered fluid coexists with the bcc crystalline structure if A 3 ^ 0. In the bcc 
arrangement of colloids the wave vectors k. J b form edges of a regular tetrahedron, and for our 
particular case a 3 = 4A 3 /(3v^6) and a 4 = 5v4 4 /8 [66]. The above results and Eqs.(f3*5J) and 
(|41|) enable us to obtain the explicit expression for the transition line in MF 

24cos(2fc b ) , s 

T = -1.063 -s w 13.1s w l.UZn p . (47) 

K 

The transition line is given by the above equation only for n p = 0(1/Z), because our 
asymptotic analysis here is restricted to s = 0(Z°). 

In the colloid limit it is more convenient to use a + = 2ct + _ as a length unit, and a^ 1 as 
a wavelength unit. To avoid confusion, the wave numbers in a^} units will be denoted by 
q. In real space the field (J4*3~j) in the case of the bcc structure is for r = (x, y, z) in suitably 
chosen coordinate frame given by 

Ap+(x, y, z) oc ip(x, y, z) oc cos ( ^7^ — ) + cos ~ ) + cos ( '^/f ~~ ) + ^ 

fQb(x-y)\ fq b (x-z)\ (q b (y-z)\ 

cos {-^-) + cos {—^-) + cos {-^-) ■ 
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FIG. 1: The surface Ap + {x,y, z) = 0, separating the regions of enhanced and depleted density 
of particles. Ap + (x,y, z) > inside the droplets. The cubic unit cell with the lattice constant 
a = 2\/2vr/(75 is shown. 

In Fig.l the surface Ap + (x, y, z) = is shown. This surface separates the regions with the 
particle density exceeding the average value from the regions of depleted particle density. 
From (J48|) we see that the lattice constant a of the bcc structure is related to the critical 
wavenumber by a = 2\/2ir/qt, « 3.6 in a + units. The distance between nearest-neighbors 
in the bcc crystal is \/67r/g& ~ 3.12cr + . Let us consider the above transition between the 



28, 



uniform phase and the bcc crystal for typical (deionized) aqueous systems studied in Refs 
29 . 130L where Z ~ 10 3 — 10 5 and A ~ 10 3 . From (fTrjj) we find that room temperature is 



T* ~ 10 1 — 10 3 , and Eq. ()47)1 gives n p ~ 10 4 — 10 8 in the uniform phase at the 
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coexistence with the bcc crystal. Such dilutions can be recognized as 'voids' in experiments. 
Moreover, in experiments the interparticle distance in the bcc crystal was found to be ~ 3a + 
29). |30j , in agreement with the results of our theory. The above behavior is not confirmed 



28, 



by simulations 



igreement w 

QQ,Q 



where gas-liquid type separation with a critical point, rather 



12 



3,0 



the size- and/or charge 



than crystallization, has been observed. However, in Refs 
asymmetry is two-three orders of magnitude smaller than in experiments. 

The asymptotic theory described above is strictly valid only in the colloid limit (A — > oo 
first and next the asymptotic behavior for Z — > oo with n p ~ 1/Z or smaller is considered). 
Beyond the colloid limit the full set of equations have to be solved. Since there are two 
eigenmodes in the general case, both eigenvalues can vanish, lea ding to two spinodal lines 



associated with two phase transitions. As will be shown in Ref. 



671 the other spinodal is 



associated with a gas-liquid separation. Beyond MF the crystallization may be preempted 



24, 



26, 561. The latter transition 



by the gas-liquid separation, as is the case for the RPM 
corresponds to vanishingly low values of n p and T* in the limit Z, X — > oo, therefore only 
the crystallization survives in the colloid limit. Thus, although predictions of our theory 
in the colloid limit disagree with the results of simulations (obtained beyond that limit) it 
is plausible that the results of the full theory will agree with simulations for appropriate 
values of the asymmetry parameters and for the corresponding regions in the phase diagram 



(n p ,T*). It is worth noting that in snapshots shown in Ref 
are clearly distinguishable. 



12j clusters separated by 'voids' 



B. Gaussian correlation functions 



For temperatures lower than that given in Eq. (j47|) the periodic structure is stable. This 
suggests effective attractions between like-charged macro ions at the distances r ~ 3a + , at 
least near the transition to the crystalline phase. In fact already the experimental discovery 



3C | inspired a debate on the origin of the 



30 



36 



38 



3, 41, 3- 



of void-crystal coexistence and other anomalies 
effective attraction between like charged particles 

In the mesoscopic theory instead of effective interactions between the macroions in the 
uniform phase we consider the correlation function G ++ (ri — r2) defined in Eq. (jl2j) . and 
related to the pair correlation function according to Eq. ()29j) . Maxima of G ++ (ri — r2) 
indicate increased probability of finding a pair of colloid particles at the corresponding 



17 



positions. In the lowest order, Gaussian approximation (neglected VL int in Eq. (j2*U|l ) the 
correlation functions are given in Eq. (j28|) . and in the colloid limit (i.e. for n p ~ 1/Z and 
A, Z — > oo) we find 

-i -l 

An cos q 



g + M = Tz 



s + 



and 



where 



7r Sq 2 



Sq 2 



47rcos(g/2) - 



(49) 
(50) 

(51) 
(52) 



S = vrT7(24s), 

and terms 0(Z~ 2 ) have been neglected. Note that the g-dependent parts of the correlation 
functions are all of the same order 0(Z~ l ). The g-dependent parts of the functions ()49|) - 
(fBTj) multiplied by 4Z/3* are independent of Z, and depend on the thermodynamic state only 
through S. In the colloid limit the functions AZj3*G a p assume universal shapes along the 
straight lines (|52|) in the phase diagram (s,T*), at least in the Gaussian approximation. 

From the above result and from Eq.(|18J) it follows that for s = 0(Z°) the corresponding 
Sf-functions (Ea. tfZty ) are g ++ = O(Z), g + _ = 0(Z°) and = 0(Z- 1 ). Note the strong 
dependence of these functions on the charge asymmetry resulting from the difference in the 
number densities in the charge-neutral system. 

In real-space representation G a p{r) can be obtained by residue method 

HQ- All the 

functions have the same denominator, hence the same poles determine the decay lengths 
and (where applicable) the period of damped oscillations. The form of G ++ (q) is similar 
to the form of charge- density correlation function in the RPM, and the latter was studied 



in RefJ65l. From the results of RefJ65l it follows that for S > Sk ~ 11.8, where S = Sk is 
known as the Kirkwood line 3], there are two imaginary poles ia% and ia>2 in the upper 



half of the complex plane, and 

rG aP {r) 



A {l) e~ 



■air | a( 2 ) p - 



For S < Sk there are two conjugate complex poles, q\ = a\ + iao and q2 
and [65j 

rG aP {r) = A^sin^r + 6)e- a ° r . 



(53) 

-ql with ao > 0, 
(54) 
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The poles can only be found numerically, except near the spinodal line 65]. Following the 
analysis of Refl65l we find the characteristic lengths and the amplitudes. The correlation 



functions AZ(3*G a p are shown in Figs. 2-4 for three different regimes. In Fig. 2 we show the 
correlation function above the Kirkwood line, i.e. for very dilute systems or for very high 
temperatures. Fig.3 corresponds to S < Sk, i.e to denser system and/or lower temperatures, 
but far from the phase coexistence. Finally, in Fig. 4 we show the correlation functions at 
the coexistence with the bcc crystal. 

From Eq.()53J) it follows that for r ^> a^ 1 (where 02 > a\ and the numbers are well 
separated) we recover the well known Yukava-type decay of correlations (see Fig.2), expected 
for very dilute systems. In this regime the effective interactions between the like-charge ions 
are purely repulsive, as in the DLVO theory. 

In Fig.3 we see a qualitative change in the shape of the correlation functions, which in this 
part of the phase diagram exhibit oscillatory decay (|54|). Let us first analyze 4Z/3*G + _(r) 

and 4Z/3*G (r). The 4Z/3*G^ (r) assumes a maximum for r rs cr + /2 and then decreases 

rather slowly for increasing r. For 2a + < r < 3a + AZ/3*G + -(r) is negative, and assumes 

a minimum for r »s 2.5cr+. At the same time AZj3*G (r) > for 1.5 04. < r < 3cr + , 

and assumes a positive maximum for r « 1.8<r + . This means that the counterions are 
preferably separated by distances 1.5cr+ < r < 3cr + , i.e. when there is more than enough 
room for a colloid particle to be located between them. The behavior of the two correlation 
functions suggests a tendency for ordering in a structure where a diffuse cloud of counterions 
is formed around the colloid particle. The cloud of counterions extends to the distance from 
the center of the colloid r « 2a + . For the distance from the center of the colloid particle 
r > 2cr + the density of the counterions is depleted compared to p$_. Let us turn to the 
4Z/3*G ++ (r). It assumes a positive maximum at r w 3a + , indicating preferable location 
of the corresponding pair of ions at such distances, consistent with formation of the cloud 
of counterions between them. Note that the clouds surrounding the two colloid particles 
separated by r « 3<7+ overlap weakly in a small region around half the distance between the 
particles. The maximum of 4Zf3*G ++ (r) in Fig.3 is only slightly larger from zero, and the 
tendency of the colloids to be separated by such a distance is very weak. 

Let us finally analyze the correlations in the uniform phase at the coexistence with the 
bcc crystal. Note first that the phase coexistence (}4Tj) occurs quite close to the spinodal 
line (jSHJ), where the amplitudes of the correlation functions diverge, and for S — > behave 
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FIG. 2: Correlation functions 4Z/3*G aj g(r) (dimensionless) in the uniform phase for S = 15, i.e. for 
very strong dilutions. The inverse decay lengths in Ea. H53|) are a\ = 1.265 and a<i = 3.19. Distance 
is in units of the particle diameter. As discussed in sec. Ill CI results of the mesoscopic theory for 
r < a are not expected to be correct. 
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FIG. 3: Correlation functions AZ(3*G a p{r) (dimensionless) for 5 = 5, i.e. in the uniform phase at 
larger densities, but still far from the transition to the bcc crystal. The characteristic lengths in 
Ea.(|54[) are ao = 1-45 and a\ = 1.72. Distance is in units of the particle diameter. As discussed in 
sec. Ill CI results of the mesoscopic theory for r < a are not expected to be correct. 
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FIG. 4: Correlation functions AZf3*G a p{r) (dimensionless) for S = 1.712, i.e. in the uniform 
phase at the transition to the bcc crystal. The characteristic lengths in Ea.(|54|l are ao = 0.3 and 
Oi\ = 2.43. Distance is in units of the particle diameter. As discussed in sec. Ill CI results of the 
mesoscopic theory for r < a are not expected to be correct. 
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as ~ (S — Sb)^ 1 / 2 in the Gaussian approximation 65[ (S is defined in Eq. flo^|) . and Sj, = 
nT b * / (24s)) . Thus, the amplitudes of the correlation functions are large at the coexistence 
with the crystal, and the fluid phase is strongly structured. To investigate this structure 

in more detail, consider first 4Zf3*G + _(r) and 4Zf3*G (r). We find a rather large value 

of 4Z(3*G + _(r) for r pa a + /2, namely 4Zj3*G + -(o + J) pa 0.737, and a rather fast decay 
of G + -(r) for increasing r; 4Z(3*G + -(r) < for 1.2er+ < r < 2.5a + . The counterion 

correlation function 4Z(5*G (r) assumes a maximum for r pa cr + , i.e. when the counterions 

are located at the opposite sides of the colloid-particle surface. From the above observations 
we can deduce that the cloud of counterions becomes much denser and thiner, and is closely 
attached to the particles. Consider now 4Z(3*G ++ (r). Positive value of 4Zf3*G ++ (a + ) 
suggests effective attraction between the like charge macroions at the distance of the closest 
approach. Similar result was also found in Reffl2. This can only be possible if the point-like 
counterions are attached to the colloid surface, consistent with the formation of a thin and 
dense layer of counterions around each colloid particle. The subsequent, positive maximum 
of 4Z(3*G ++ (r) occurs at r pa 3cr + . This maximum is much higher than away from the 
phase transition (Fig. 3), and suggests rather strong tendency for location of colloids at such 
distances. Note that the clouds of counterions around the colloid particles separated by the 
distance r pa 3<j + do not overlap. The second maximum of 4Z(3*G + -(r) at r pa 3<7+, i.e. at 
a similar distance, shows again that each colloid particle is surrounded by a dense and thin 
cloud of counterions. 

The correlation functions g a p were obtained in simulations for Z = 10, A = 19. 
These asymmetry parameters are too small for our asymptotic regime A, Z — > oo, and the 
correction terms in Eqs. (j49j) - ([5T| may well be of the same order as the terms which in the 
asymptotic regime dominate. This may be an important source of discrepancy between our 
theory and simulations. Also, we only obtained the Gaussian correlation functions in the 
local-density approximation. The main discrepancy between our results and the results of 
simulations concerns the positions of the maxima of the correlation functions. Note, however 
that in experiments the colloidal crystals are formed in the case of very strong size and charge 
asymmetry two-three orders of magnitude larger than studied in simulations jl^ . 
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IV. SUMMARY AND DISCUSSION 



We have developed a mesoscopic theory for the PM with arbitrary size- and charge 
asymmetry. Our theory allows for systematic studies of phase transitions and structure for 
any charge- and size ratio. The results obtained in MF and Gaussian approximations can 
be improved by adding fluctuation-corrections obtained in perturbation theory. 

Explicit results for phase transitions and structure in the uniform phase were obtained in 
the colloid limit (A — > oo first, and next the asymptotic regime of Z — > oo with n p = 0(1/ Z) 
is considered) in the MF approximation. We found a coexistence of a very dilute phase 
with the bcc crystal formed by the colloid particles. The lattice constant was found to be 
a ~ 3.6cr + . Very strong dilution of colloids in the 'gas' phase, structure of the crystalline 
phase and the lattice constant agree with experimental results. The correlation functions 
Gap( T ) f° r density deviations of the species a, (3 — ± at the distance r show the known 
monotonic decay for large values of S defined in Eq. fl52|) (high temperatures T* and/or low 
densities s). For decreasing 5* the short-range order in the uniform phase increases. 

In this work the analysis of the colloid limit is restricted to the MF approximation. 
Inclusion of fluctuations will certainly change the quantitative results, in particular the 
location of the phase transition. We expect that the fluctuations do not play a dominant 
role in the colloid limit, but the role of fluctuations certainly deserves attention in future 
works. In the full theory two spinodal lines occur, and coupling between the fields p* + and 
p*_ in AQ MF may lead to an increased role of fluctuations. By analogy with the RPM 



we expect that for not too large values of A, Z and/or for volume fractions larger 
than ~ 1/Z, fluctuations may induce significant shifts of the spinodal lines, including the 
change of metastable transitions into stable ones and vice versa. Hence, beyond MF the 
crystallization may be preempted at low concentrations and temperatures by the gas-liquid 
type separation for certain values of Z and A. The role of fluctuations for different asymmetry 
parameters will be studied in future works. 

We should emphasize that the foundations of the mesoscopic description and the asymp- 
totic analysis for large asymmetry are based on first-principle considerations rather than 
having been fit to the results of experiments. Mesoscopic field theories turned out to be 
appropriate for a description of a weak ordering, including a weak crystallization. Because 
the unit cell of the experimentaly observed bcc crystal |2Ja] is rather large, one may expect 
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that the corresponding transition is not associated with close packing. The nearest-neighbor 
distance in the ordered structure corresponds to the second maximum in the corresponding 
correlation function in the uniform phase close to the phase coexistence. As the results of 



54, 



651 ] . it is plausible 



mesoscopic theories are quite accurate down to such distances 
that in this particular case our theory yields correct results on a semiquantitative level. 
Obviously, our mesoscopic field theory has its limitations, and the structure for distances 
w a cannot be correctly reproduced, as is also the case in the commonly accepted Landau- 
Ginzburg- Wilson and Brazovskii theories. Our theory should be considered as a contribution 
to the discussion concerning the thermodynamics and structure in the charged colloidal sys- 
tems. Both the experiments and our theory show the formation of the bcc structure with 
a large unit cell. To confirm that this is a real phenomenon it is desirable that microscopic 
theories and/or simulations yield similar results. 

Acknowledgments 

The work of AC and WTG was partially funded by the KBN grant No. 1 P03B 033 26. 
We also gratefully acknowledge the support of the Division of Chemical Sciences, Office of 
the Basic Energy Sciences, Office of Energy Research, US Department of Energy. 



[1] G. Stell, K. Wu, and B. Larsen, Phys. Rev. Lett. 37, 1369 (1976). 

[2] M. E. Fisher, J. Stat. Phys. 75, 1 (1994). 

[3] C. W. Outhwaite, Cond.-Matter Phys. 7, 719 (2004). 

[4] G. Stell, J. Stat. Phys. 78, 197 (1995). 

[5] L. Onsager, J. Am. Chem.Soc. 86, 3421 (1968). 

[6] G. Stell and J. Lebowitz, J. Chem. Phys. 49, 3706 (1968), see esp. Eqs. (2.9), (2.12) and 
(2.13). 

[7] G. Stell, , sec. 2.2 in Ref. 4. 

[8] J. M. Romero-Enrique, G. Orkoulas, A. Panagiotopoulos, and M. Fisher, Phys. Rev. Lett. 85, 
4558 (2000). 

[9] Q. Yan and J. de Pablo, Phys. Rev. Lett. 86, 2054 (2001). 



25 



[10] Q. Yan and J. de Pablo, Phys. Rev. Lett. 88, 95504 (2002). 

[11] A. Panagiotopoulos and M. Fisher, Phys. Rev. Lett. 88, 45701 (2002). 

[12] D. Cheong and A. Panagiotopoulos, J. Chem. Phys. 119, 8526 (2003). 

[13] J. Rescic and P. Linse, J. Chem. Phys. 114, 10131 (2001). 

[14] P. Linse, Philos. Trans. R. Soc. London, Ser.A 359, 853 (2001). 

[15] Y. V. Kalyuzhnyi, M. F. Holovko, and V. Vlachy, J. Stat. Phys. 100, 243 (2000). 

[16] D. M. Zukerman, M. E. Fisher, and S. Bekiranov, Phys. Rev. E 64, 11206 (2001). 

[17] M. N. Artyomov, V. Kobelev, and A. B. Kolomeisky, J. Chem. Phys. 118, 6394 (2003). 

[18] J.-N. Aqua and M. E. Fisher, Phys. Rev. Lett. 92, 135702 (2004). 

[19] F. Raineri, J. Routh, and G. Stell, J. Phys. IV (France) 10, 99 (2000). 

[20] G. Stell, in New Approaches to Problems in Liquid-State Theory, edited by C. Caccamo, J.-P. 

Hansen, and G. Stell (Kluwer Academic Publishers, Dordrecht, 1999). 
[21] G. Stell, Phys. Rev. A 45, 7628 (1992). 
[22] G. Stell, , sec.2.1.2 in Ref. 4. 

[23] A. Ciach and G. Stell, J. Mol. Liq. 87, 253 (2000). 

[24] A. Ciach and G. Stell, Physica A 306, 220 (2002). 

[25] A. Ciach and G. Stell, Int.J. Mod. Phys. B 21, 3309 (2005). 

[26] A. Ciach and G. Stell, Phys. Rev. E 70, 16114 (2004). 

[27] A. Ciach, Phys. Rev. E 70, 046103 (2004). 

[28] N. Ise et al, J. Chem. Phys. 78, 536 (1983). 

[29] N. Ise, T. Konish, and B. Tata, Langmuir 15, 4176 (1999). 

[30] A. K. Arora and B. V. R. Tata, Adv. Colloid Interface Sci. 78, 49 (1998). 

[31] M. Leunissen et al, Nature 437, 235 (2005). 

[32] B. V. R. Tata, E. Yamahara, P. V. Rajamani, and E. Ise, Phys. Rev. Lett. 78, 2660 (1997). 

[33] K. Ito, H. Yoshida, and N. Ise, Science 263, 66 (1994). 

[34] B. Derjaguin and L. D. Landau, Acta Physicochim. URSS 14, 633 (1941). 

[35] E. J. W. Verwey and J. T. G. Overbeek, Theory of the Stability of Lyophobic Collids (Elsevier, 

Amsterdam, 1948). 
[36] P. B. Warren, J. Chem. Phys. 112, 4683 (2000). 

[37] M. C. Barbosa, M. Deserno, C. Holm, and R. Messina, Phys. Rev. E 69, 51401 (2004). 
[38] L. Belloni and O. Spalla, J. Chem. Phys. 107, 465 (1997). 

26 



[39] R. D. Groot, J. Chem. Phys. 94, 5083 (1991). 

[40] R. Messina, C. Holm, and K. Kremer, Phys. Rev. Lett. 85, 872 (2000). 

[41] E. Allahyarov, I. D. Amico, and H. Lowen, Phys. Rev. Lett. 81, 1334 (1998). 

[42] H. Lowen, P. A. Madden, and J.-R Hansen, Phys. Rev. Lett. 68, 1081 (1992). 

[43] Y.-X. Yu, J. Wu, and G.-H. Gao, J. Chem. Phys. 120, 7223 (2004). 

[44] I. S. Sogami and N. Ise, J. Chem. Phys. 81, 6320 (1984). 

[45] R. van Roij, M. Dijkstra, and J.-P. Hansen, Phys. Rev. E 59, 2010 (1999). 

[46] M. N. Tamashiro and H. Schiessel, J. Chem. Phys. 119, 1855 (2003). 

[47] R. R. Netz and H. Orland, Europhys. Lett. 45, 726 (1999). 

[48] R. R. Netz and H. Orland, Eur. Phys.J. E 1, 67 (2000). 

[49] J. M. Caillol, J. Stat. Phys. 115, 1461 (2004). 

[50] A. Ciach and G. Stell, J. Chem. Phys. 114, 382 (2001). 

[51] J. Jiang et al, J. Chem. Phys. 116, 7977 (2002). 

[52] S. A. Brazovskii, Sov. Phys. JE TP 41, 8 (1975). 

[53] R. Evans, R. L. de Carvalho, J. R. Henderson, and D. C. Hoyle, J. Chem. Phys. 100, 591 
(1994). 

[54] R. L. de Carvalho and R. Evans, Mol. Phys. 83, 619 (1994). 

[55] G. Stell, in New Approaches to Problems in Liquid-State Theory, edited by C. Caccamo, J.-P. 

Hansen, and G. Stell (Kluwer Academic Publishers, Dordrecht, 1999). 

[56] A. Ciach and G. Stell, Phys. Rev. Lett. 91, 60601 (2003). 

[57] Y. Shim, M. Y. Choi, and H. J. Kim, J. Chem. Phys. 122, 044510 (2005). 

[58] A. Z. Panagiotopoulos and S. Kumar, Phys. Rev. Lett. 83, 2981 (1999). 

[59] A.Diehl and A.Z.Panagiotopoulos, J. Chem. Phys. 118, 4993 (2003). 

[60] A. Diehl and A. Z. Panagiotopolous, Phys. Rev. E71, 046118 (2005). 

[61] C. V. F.Bresme and J.L.F.Abascal, Phys. Rev. Lett. 85, 3217 (2000). 

[62] C. Vega, J. Abascal, C. McBride, and F. Bresme, J. Chem. Phys. 119, 964 (2003). 

[63] A. Ciach and G. Stell, J. Chem. Phys. 114, 3617 (2001). 

[64] J. L. Lebowitz and J. S. Rowlinson, J. Chem. Phys. 41, 133 (1964). 

[65] A. Ciach, W. T. Gozdz, and R.Evans, J. Chem. Phys. 118, 3702 (2003). 

[66] L. Leibler, Macromolecules 13, 1602 (1980). 

[67] A. Ciach, W. T. Gozdz, and G. Stell, to be published . 



27 



[68] L. Blum, Mol. Phys. 30, 1529 (1975). 

[69] L. Blum and J. S. Hoye, J. Phys. Chem. 81, 1311 (1977). 

[70] K. Hiroike, Mol. Phys. 33, 1195 (1977). 

[71] I. Yukhnovskii, Sov.Phys. JETP 34, 263 (1958). 

[72] I. R. Yukhnovskii, Phase Transitions of the Second Order, Collective Variable Methods (World 

Scientific, Singapore, 1978). 

[73] O. Patsahan, Cond.Mat.Phys. 7, 35 (2004). 

[74] J. Hoye and G. Stell, J. Stat. Phys. 89, 177 (1997). 



28 



